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Abstract In the ultra-cold regime Van der Waals forces between neutral atoms can be rep- 
resented by short range effective interactions. We show that universal low energy scaling 
features of the underlying vdW long range force stemming from two photon exchange im- 
pose restrictions on an Effective Field Theory without explicit photons. The role of naively 
redundant operators, relevant to the definition of three body forces, is also analyzed. 
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1 Introduction 

From a fundamental QED point of view the underlying mechanism responsible for van 
der Waals (vdW) forces corresponds to two photon-exchange (see e.g. 1 1 ] and references 
therein). As compared to the short range and exponentially suppressed chemical bonding 
forces on sizes about a few Bohr radii, vdW forces are long range. For interatomic separa- 
tions flB ^ r ^ he/ AE, two photons are exchanged in a short time ~ 2r/c while transitions 
with excitation energy AE take a much larger time ^ Ih/ AE, yielding the potential 

^W = -76 -7^-710 W 

where C„ are the dispersion coefficients which are accurately known for many diatomic 
systems (see e.g. a compilation in (2|)- The vdW length R = (MCe/Zj characterizes the 
size of the forces. For such potentials, low energy scattering with kR '\& dominated by 
S-waves which phase-shift, 5(,{k) , fulfills the effective range expansion (ERE) [SJ 

kcoi5Q{k) = -— + \ rok^ + V2k'^ log{k^R^ ) + ... (2) 
«() 2 

where ocq is the scattering length, and ro is the effective range. Note that for this potential 
the long-range character stars at ff{k'^) due to the logarithmic piece. 
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2 Low energy Scaling of vdW forces 

Remarkably, the effective range was computed analytically HJSj when C„>8 = yielding 

ro 16r(5/4)^ 4R 4r(3/4fR^ R R^ 

— = ; ' ' --— + 7 = 1.395- 1.333— +0.6373^ (3) 

R 371 3 Oo 3k Oq 

The scaling of the effective range ro in the vdW length and the quadratic l/oo behaviour is 
just a particular case of a more general result |6| (see also |7| in these proceedings). In fact, 
the dominance of the leading long distance Cg term tacitly assumed in Refs. |4l5| was to be 
expected a priori by suitably re-writing higher order Cg, Cm contributions in vdW units 

R^MV{r) = -{R/rf [l + gi{R/rf + g2{R/ r)U . . .] (4) 

where gi ~ 10^^ and g2 ~ 10^"* for many homonuclear diatomic systems. Thus one expects 
that even for kR ^ I higher order corrections are negligible despite the strong divergence at 
short distance. These expectations are indeed met a posteriori on the light of about a hun- 
dred calculations based on phenomenological potentials (6]. This result not only favours the 
view that these rather simple approaches based on the leading vdW forces are phenomeno- 
logically sound but also shows that a huge reduction of parameters takes place suggesting 
that atoms in the ultra-cold regime can indeed be handled without much explicit reference to 
the underlying electronic structure of atoms. The scaling universal relation, Eq. ([3]l, allows 
for a quite general discussion on effective interactions in vdW units, as we advance here. 



3 Effective Short Distance Potentials 

In the ultra-cold regime, i.e. for extremely long de Broglie wavelengths much larger than the 
vdW scale, X = \ /k^ R , one expects the long range character to become largely irrelevant, 
keeping the first two terms in Eq. l|2]l. Thus, one might want to represent the vdW potential 
by an effective potential with a finite range, rc, featuring the truncated ERE, Eq. and 
dismissing any explicit reference to the underlying photon exchange. However, even at very 
low energies, causality arguments provide the shortest possible value for rc, which for vdW 
forces yields rc > 0.6/? |2|. Using for illustration a square well (SW) potential with range rc 
and depth Vb, V'eff(^) = — Vb0(''c — r), one obtains 

tsn y/M%rc r 1 

0^0 = ''c — : ''0 = rc 



1- 

OarcMVi) 3al 



(5) 



Reproducing Eq. (O is not possible for a common potential. Indeed, the sign of the l/a^ 
term is just opposite, so that for small ocq we cannot represent the interaction by this short 
range potential. On the contrary, for large scattering lengths ao 3> ^ we obtain = 1.395^ 
and Vo = 7r^/(4r^M). In terms of volume integrals, 

Cq = j d\v,n{iL) , C2 = -^-j d\,^V,n{iL) (6) 

one gets MCf^ /R = -14.41 and MCl^ /R^ = 2.80. If we use instead a delta-shell (DS) 
potential Veff('') = -Vorc5{r- rc) we get for Oo > i? the results MCq^/R = -13.15 and 
MC2^ /R^ = 2.40, not far from the SW estimate. This suggests using a formulation based 
directly on the constants Co and €2- Note that while a C4 exists for these short distance 
potentials, the original vdW potential yields a divergence, in harmony with the observation 
that the ERE for short range potentials differs at (?{k'^) from the vdW expression, Eq. ([3]l. 
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4 Effective Field Theory 

The EFT approach has often been invoked to highhght universal features of ultra-cold few 
atoms systems (for reviews see e.g. 1 8 , 9 1). We re-analyze it on the light of the universal and 
extremely successful scaling relation, Eq. ([3]l. For definiteness, we consider the Galilean 
invariant Lagrangian density 1101 expanded in composite Bosonic spinless field operators 
with increasing energy dimensions and including multi-particle interactions, 

Co ,.,,,....,2 C2 .,^,,2 ,.,^,3 



^ = + 2^)'''- -f^'f 'fr - y [v(v^V)]' - -^(v^V)' + • • ■ (7) 

Here Co, C2 and Dq are low energy constants which are fixed from few body dynamics. 
Using Feynman rules in the two-body sector one derives a scale dependent and momentum 
truncated self-adjoint pseudo-potential in the CM system (k and k' are relative momenta) 

(k'|y|k)= [Co + C2{k^+k'^) + ...]e{A-k)e{A-k'). (8) 

The cut-off A is introduced here to handle the power divergent integrals arising in the scat- 
tering problem, which in terms of the Lippmann-Schwinger (LS) equation becomes 

implementing unitarity for p < A. Using the potential of Eq. ^ the LS Eq. (|9]l reduces to a 
system of algebraic equations which solution is well known (see e.g. Ref. 111,1 ) yielding 

1 4(-2c^+90;r4-H5(3co+2c2);r^) 



oqA 9k[cI-IOcoK 

_ 16{4 + 12k^C2+971'^) 12C2{C2 + 1271^) 1 3c27l {c2 + 12k^) 1 

~ n{c2 + 6n^f {c2 + 67z^f «oA (^c2+67:^f «oA 



where cq = MACo, C2 = MA^C2. By eliminating Co in terms of ao we have written ro in a 
form similar to Eq. ([3]l. This leads for any cut-off A to the mapping {(XQ,rQ} (Co,C2). For 
C2 = one gets ro = 4/7iA which for ao > ^ yields AR = 0.91 and MCq/R = -21.6 from 
matching the scattering length and the effective ranges r^f^ = ''o'^- The cut-off dependence 
for C2 ^ can be looked up at Fig. [T]in vdW units and for the specific case ocq/R = 10 
where a weakly bound state takes place. As we see there is a clear stability plateau in the 
region A ~ k/{2R) illustrating the basic point of the EFT; low energy physics is cut-off 
independent within a given cut-off window which does not resolve length scales shorter 
than the vdW scale. Numerically we get MCq/R ~ -15 and MC2/R^ ~ 2 for A ~ n/{2R), 
in agreement with the previous SW and DS analysis. The values of A where the EFT low 
energy parameters diverge correspond to an upper bound above which Co and C2 become 
complex, violating the self-adjointness of the potential O and the Lagrangian, J^{x) ^ 
J^\x). Thus, off-shell two-body unitarity and hence three-body unitarity are jeopardized 
for AR > 4, despite the phase shift being real and on-shell unitarity being fulfilled. 

Direct inspection shows that a perfect matching between the vdW and the EFT effec- 
tive ranges, Eq. ^ and Eq. l llOt for any ocq is not possible. So we try out including re- 
dundant operators which are usually discarded IIIOI but are needed to guarantee off-shell 
renormalizability of the LS equation 1121 . A Galilean invariant term of the form A^ = 
— ^C2i^f^ ^f) (^td, + /2m'J is formally redundant since it can be eliminated by a 
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Fig. 1 Cut-off dependence of the EFT coefficients MCo/R when C2 = (dashed red) and MCo/P? and 
MC2IP? after Eqs. (lO) (full blue) for the case ao/R= 10. R is the vdW scale defined as i? = (MCe/S^) ' . 

field transformation i// — )• jC'',\\r{\\r^\\f) which generates additional three body forces as 
well. The new term adds a correction AV = C'2(2p^ — — k'^) /2 to the potential, Eq. ([8]l, 
vanishing on-shell. Solving the LS equation and eliminating Co in terms of Oq yields 

16 ((4 - 2c2)2 + 36;r4 +6(8c2 - 4);r2) 
roA = = 

7z{-2c2 + c'2-nn^) 
12((4-2c2)2+48c2;r2) i 37t {{c'^-lcif +4Sc27t^) 1 
(-2c2 + c2-12;r2)^ (_2c2 +4 - 12;r2)^ ^^«o ' 

where c'2 = MA^Cj appears through the combination C'2 — IC2 which cannot be completely 
eliminated by making C2 — s- C2 + 5C2. Note the accidental correlation —Ajn between the 
second and the third coefficients holding regardless on the particular regularization method. 
Perfect matching can only be achieved with complex coefficients. Minimizing the difference 
between rj^^ and r^f^ provides a reasonable range AR = 1.6— 1.8 ~ 7t/2. As we can see, 
universal two-body scaling features encoded in Eq. (O and exhibiting the underlying vdW 
(two photon exchange) nature of interactions impose severe restrictions on the EFT solu- 
tion with no explicit photonic degrees of freedom and distinguish between naively unitarily 
equivalent Hamiltonians mixing different particle number (see e.g. Ref. 1131 ). Therefore, 
these limitations are expected to play a role in the EFT analysis of three-body forces. 
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